Nonclassical characteristic functions for highly sensitive measurements 
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Abstract 

Characteristic functions are shown to be useful for highly sensitive measurements. Redistributions of motional Fock states 
of a trapped atom can be directly monitored via the most fragile nonclassical part of the characteristic function. The method 
can also be used for decoherence measurements in optical quantum-information systems. 

PACS numbers: 03.65.Bz, 42.50.Ar, 42.50.Dv 



I. INTRODUCTION 

The experimental demonstrations of photon anti- 
bunching sub-Poissonian photon statistics Q and 
quadrature squeezing Q also led to an increasing inter- 
est in practical applications of nonclassical states. An 
early example is the proposal to use squeezed light for 
enhancing the sensitivity of intcrfcromctric gravitational- 
waves detection 0] . Experiments have demonstrated the 
usefulness of squeezed light for improving interferometric 
measurements [1, 0] and spectroscopy 0]. 

Two decades after the first experimental demonstra- 
tions of the potential usefulness of nonclassical states the 
latter still play a minor role in practical measurements. 
There may be several reasons for this fact. First, the ex- 
perimental effort for generating the needed nonclassical 
states is rather high. Second, some applications, e.g. the 
use of squeezed light for optimizing the laser power in 
gravitational-wave detection, can be replaced with devel- 
opments of laser sources. Third, nonclassical states are 
usually highly fragile against losses which may substan- 
tially limit their advantages in some applications. 

The use of nonclassical states is frequently considered 
in the context of the reduction of the quantum noise in a 
certain observable below an ultimate classical noise limit. 
Examples are the use of sub-Poissonian and squeezed 
light fields for reducing the noise in direct and homodync 
photodctection, respectively. This requires to link the 
measurement principle with the observable whose quan- 
tum noise is reduced. Below we will reconsider the appli- 
cation of nonclassical states from a much broader point of 
view. When speaking about nonclassical states in the fol- 
lowing, we will only consider quantum states of a single- 
mode harmonic oscillator whose Glauber-Sudarshan P- 
function is not a probability density Q. 

The nonclassicality of quantum states can be com- 
pletely characterized in terms of measurable character- 
istic functions of phase-dependent quadratures, 
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fi (fit) being the bosonic annihilation (creation) operator 
and if is the phase parameter. To be more specific, a 
hierarchy of necessary and sufficient conditions has been 
derived that completely characterizes the nonclassical- 



ity of a given quantum state in terms of the quadrature 
characteristic function G(k, ip) [9[. A broad class of non- 
classical states can be well characterized by the rather 
simple condition of first-order nonclassicality, 

\G(k, (p)\ > Ggr(k), (2) 

stating that the absolute value of the characteristic func- 
tion exceeds, for some arguments, the corresponding 
value of the ground (or vacuum) state [l(| • 

The signatures of first-order nonclassicality are more 
general than the quantum-noise reduction of a chosen 
observable below some classical limit. The condition also 
includes features like quantum interference pd| and sub- 
Planck structures in phase space [12]. Note that the 
nonclassical effects of first and second order have been 
experimentally demonstrated for radiation fields [HI, [3] . 
The needed characteristic functions can also be observed 
for the quantized center-of-mass motion of a trapped 
ion [IE Eft 

In this paper we propose a new measurement princi- 
ple, where the quadrature characteristic function serves 
as a highly sensitive probe. It makes use of the fact that 
the nonclassical signatures of the quadrature character- 
istic function are more fragile with respect to dissipation 
than other nonclassical effects, such as sub-Poissonian 
statistics and squeezing. This provides a new tool for the 
highly sensitive diagnostics of decoherence effects, which 
is of great interest for quantum information processing. 

The paper is organized as follows. In Sec. II we con- 
sider the detection of the characteristic function of the 
quadrature distribution for the motion of trapped ions 
and for propagating radiation fields. The decoherence is 
caused in both cases, for example, by a thermal reservoir. 
Section III is devoted to the use of the most pronounced 
nonclassical features of the characteristic functions for 
highly sensitive measurements. A brief summary is given 
in Sec. IV. 



II. DECOHERENCE IN TERMS OF CHARAC- 
TERISTIC FUNCTIONS 

Let us consider a nonclassical state p(0) of a bosonic 
mode prepared at the initial time t = 0. Its further evo- 
lution is caused by the dissipation to be analyzed, leading 



1 



to the state p(t). Eventually, the quadrature character- 
istic function G(k,t,ip) is measured. 

In the case of a trapped ion the measurement of 
G(k,t,ip) is performed as follows. An electronic transi- 
tion is driven simultaneously on the red and the blue mo- 
tional sidebands in the resolved sideband regime, which 
is described by the interaction Hamiltonian 



Tr {pexp(^d t — £*a)}, has the solution [ljj| 

Xfot) = exp{-(n+l/2)|£| 2 [l-exp(-2 7 i)]} 

x x (£exp(- 7 t),0), (7 ) 

where x(£., 0) is the Wigner characteristic function of the 
initial quantum state. From this result, the observable 
quadrature characteristic function, 
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where A i:j = \i)(j\ = 1, 2) is the electronic flip oper- is easily derived - 



ator and f2 is the effective Rabi frequency. Most impor- 
tantly, it is proportional to the quadrature operator x v of 
the center-of-mass motion, the phase ip being controlled 
by the phase difference of the driving lasers. The total 
state of the ion is g(t) = p(t) (g> a(t), with a properly pre- 
pared electronic state a(t). At time t the interaction Q 
is switched on for the interaction time r. The observation 
of the occupation an(t + r, tp) of the electronic ground 
state directly yields the characteristic function G(k, t, tp) 
of the quadrature distribution (l5j : 
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where the interaction ^ leads to the scaling k = 2|0|t. 

The incoherent and coherent occupations, a^ 1 "^ and 

a[ c ° h ^ , are measured with the electronic preparations 
CTn(t) = 1 and crn(t) = |cr 12 (£) | = h, respectively. The 
electronic-state occupations in Eq. g]) are detected with 
almost perfect efficiency, by testing (at time t+r) a tran- 
sition from the state |1) to an auxiliary state for the ap- 
pearance of fluorescence [l6| . 

For a radiation field the characteristic function can be 
sampled in balanced homodyne detection [l3|, [l4| , 
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Here x v> j(t), j = 1, . . . , N v , is the set of data recorded in 
balanced homodyning for each setting of (p. The quadra- 
tures describe now the radiation mode at the time t. 
Thus they carry the information on the dissipative in- 
teraction contained in the radiation-state p(t). 

To illustrate the idea, we deal with a simple model of 
decoherence caused by a thermal bath of mean occupa- 
tion number h. The density operator p in the interaction 
picture obeys the master equation 



~rP = 7 ( n + 1) [2apa^ — a)ap — pa^a] 
at 

+7n [2a) pa — aa) p — paaJ] , 



(0) 



with 7 being the damping rate. The resulting equa- 
tion for the Wigner characteristic function, t) = 



III. HIGHLY SENSITIVE DETECTION 
A. Decoherence of a Trapped Ion 

Let us first consider a trapped atom which is initially 
in the number state |m), which can be realized in exper- 
iments fl7j . The motional-state redistributions caused 
by the reservoir lead to strong modifications of the non- 
classical signatures of the characteristic function, which 
represents our highly sensitive probe. Note that the ob- 
served decoherence of a Raman-driven trapped ion [l7j 
is not completely understood yet. Although dephasing 
mechanisms could be identified [3], a deeper insight in 
the role of motional states is still required. 

The Wigner characteristic function of the number state 
\m) is given by 



Xm(6=£m(|?| 2 )exp(-|£|72), 



(9) 



where L m (x) is a Laguerrc polynomial of order m. 
Since the state |m) is phase independent, we may write 
G m (k,ip) = G al (k). In Fig. 1 wc show the quadrature 
characteristic functions, G m (k) = G m (k,t = 0), as func- 
tions of k, for m = 9,10,11. Clearly, the first-order 
nonclassicality condition © is fulfilled for all the shown 
number states. For our examples, the first-order non- 
classical effect is most pronounced at the outermost local 
cxtrcmum, where |G m (fc)| — G gr (fc) becomes maximal. 

Now we consider the time evolution of the character- 
istic function, G m (fc,t), caused by the thermal reservoir. 
The initial preparation of a Fock state |m) allows one 
to distinguish motional-state redistribution effects from 
dephasing effects. Using Eqs ([7]) to we obtain 



G m (k,t) 



exp {-n[l - cxp(-27i)]fc 2 } 

xL m (k 2 e- 2lt ) exp (-fc 2 /2). (10) 



This function can be monitored for a chosen time t, as 
discussed in connection with Eq. 

To get more insight into the time evolution, we consider 
the time derivative of the characteristic function at t = 0, 



G m (fc,0) = 2 7 [lW# 2 ) - nL m (k 2 )} k 2 e- k2 ' 2 

It strongly depends on both k and n, for m = 
maximum of the time derivative occurs for k = km 



(11) 
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FIG. 1: Characteristic functions G m (fc) versus k for number 
states m — 10 (full line), 9 (dashed-dotted) and 11 (dotted), 
together with the classical limit G SI (k) = e~ k ^ (dashed). 

Hence this argument of the characteristic function will be 
of particular interest to detect the decoherence in a sensi- 
tive way. Comparing with Fig. 1 , the most sensitive reac- 
tion on the motional-state redistributions occurs around 
the outermost maximum of the initial characteristic func- 
tion, where the nonclassical features arc dominant. This 
reflects the expected high fragility of nonclassical effects. 

For the used Fock state it is interesting to compare the 
time evolution of the characteristic function with that of 
the sub-Poissonian statistics. For the Mandel parameter, 
Q = (((An) 2 ) — (n))/(n), which measures the deviation 
from the Poissonian statistics, we get 

_ (ft 2 - 2nm - m)e~ 47 * + 2n(m - n)e~ 27 * + n 2 
QW = me -2 7 t + n(i_ e -27t) • 

(12) 

Negative values of Q characterize a nonclassical quantum 
state showing sub-Poissonian statistics. 




FIG. 2: The time evolution is shown for the occupation prob- 
ability pio(£) of the state \m = 10), the normalized mean 
excitation number n(t), the Mandel parameter Q(t), and the 
normalized characteristic function gio(fcmax, t), for n = 1. 

In Fig. 2 we show the time evolution of the normalized 
characteristic function, gio(k,t) — Gio(fc,i)/Gio(£;,0), 
for k ~ fc max and n = 1. For comparison, the time 



evolutions are also shown for the Mandel Q parame- 
ter, the occupation probability pio(t) of the initially pre- 
pared Fock state \m = 10), and the normalized mean 
motional-state excitation n(t), n(t) = (h(t)) / (n(0)} with 
(h(t)) = me~ 27 * + h(l — e 27 *). It is clearly seen that 
the nonclassical characteristic function shows the fastest 
decay and evolves much faster than the nonclassical prop- 
erty described by the Mandel Q parameter. In fact, the 
nonclassical part of the characteristic function decays 
even faster than the occupation probability of the ini- 
tially prepared Fock state. Thus it yields a highly sensi- 
tive means of detection, since even tiny dissipation leads 
to a noticeable effect on that part of the characteristic 
function. 

How can we explain the highly sensitive behavior of 
the characteristic function in its outermost maximum? 
The quadrature characteristic function, 

oo 

G 10 (k,t) = J2Pn(t)G n (k), (13) 

n=0 

can be written as a sum of characteristic functions of 
number states weighted with the occupation probabili- 
ties p n (t). From the master equation it follows that 
at the very beginning of the time evolution the only non- 
vanishing occupation probabilities are pg(t), pn(t) and 
Pw{t) = 1 — P9(t) — Pn{t)- It is seen from Fig. 1 that 
Gg(fc) and Gu(fc) attain for k = fc max roughly the same 
absolute value as Gio(fc), but with opposite sign. Thus, 
according to Eq. (fl3|) . the increase of pg(t) and pu{t) 
leads to a faster decay of Gio(k, t) compared with pio(t). 
The decay of Gio(k max ,t) is roughly twice as fast as the 
decay of Pw(t). Hence the detection of the characteris- 
tic function is both simpler and more sensitive than a 
number-state measurement, even though the latter is di- 
rectly related to the prepared number state. 

For comparison we may look at the characteristic func- 
tion for other fc-values, which are within the classical re- 
gion (such as k f=a 1) or which violate the classical limit 
only slightly (e.g. k « 2). Then the temporal evolution is 
significantly slower. For k « 1 the neighboring functions 
are almost equal: Gg(fc) ~ Gio(fc) ~ Gn(fc), cf. Fig. 1. 
Thus increasing values of pg and pxi nearly compensate 
the decay of pio, resulting in a slow decay of G\g{k,t). 

B. Decoherence of Radiation Fields 

In the following we will consider a typical applica- 
tion of our method that could be realized for a radiation 
field. The aim is to identify small decoherence effects dur- 
ing the transmission of nonclassical light through media. 
This is of great importance for applications such as in 
quantum communication, where the transmission can be 
performed via optical fibers or through the atmosphere. 

The preparation of photon number states, at least for 
larger photon numbers, is much more difficult to realize 
as the Fock-state preparation in the ion trap. However, 
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FIG. 3: Time evolution of the normalized mean photon num- 
ber n(t), the quadrature variance s(t), and the characteristic 
function <? S v(fcmax, t) of an initial squeezed vacuum state with 
r = 1.32, for n = 1. 

we will see that squeezed light can be used as the non- 
classical radiation source for our measurement principle 
as well. Let us consider a squeezed vacuum, 

|sv) = cxp[-(r/2)a t2 + (r/2)a 2 ]|0), (14) 

with r > 0. The action of the medium is modeled by the 
thermal reservoir as before. 

Due to the contact with the reservoir, the minimum of 
the quadrature variance (ip = in Eq. ([T])) behaves like 

((Ai(i)) 2 ) min = (1 + 2n)(l - e" 2 ^) + e^e" 2 ^ 4 , (15) 

where t represents the propagation time through the 
medium. For the chosen phase, ip — 0, the nonclassical 
effect is most pronounced. The characteristic function 
for this phase is simply given by 

G SV (M) = e- fe2 « A ^» 2 >— 12 ■ (16) 



Its time derivative yields the fc-value with the maximum 
sensitivity of our method to be 

k max = V2^, (17) 

which is independent of the value of n. 

In Fig. 3 we show the time evolution of the normal- 
ized mean photon number n(t), the normally-ordered 
quadrature variance s(i), and the characteristic function 
g sv {k,t), where s(t) = (: (Ax{t)) 2 :) min /(: (Ai(0)) 2 :) min 
and g sv (k,t) = G sv (k,t)/G sv (k,0). As expected, the 
fastest decay is observed for the characteristic function at 
k = A; max . To detect this highly sensitive reaction on the 
reservoir effects, the characteristic function can be sam- 
pled, at the end of the transmission channel, via balanced 
homodyne detection. For the sensitive detection method 
under study, we only need to consider the behavior for 
k = femax an( l f° r the phase with the minimal quadrature 
variance. 

IV. SUMMARY 

In conclusion, we have shown that highly sensi- 
tive measurements can be performed by detecting the 
strongly nonclassical part of the quadrature characteris- 
tic function. The method makes use of the high fragility 
of the nonclassical effects. For example, the direct ob- 
servation of the characteristic function can monitor the 
redistribution of the motional-state occupations of an ini- 
tially prepared Fock state of a trapped ion. Sensitive op- 
tical decoherence measurements can be realized by using 
squeezed light. The method may be useful for highly sen- 
sitive noise control in quantum information systems. It is 
based on a universal measurement principle, that can be 
used with different types of initially prepared nonclassical 
states. 
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